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Abstract: 

Dispersion relations and amplitudes of collective pionic modes are derived in a vr + 
NN~'^ + AA^~^ model for use in transport descriptions by means of a local density 
approximation. It is discussed how pionic modes can be converted to real particles 
when penetrating the nuclear surface and how earlier treatments can be improved. 
When the surface is stationary only free pions emerge. The time- dependent situation 
is also addressed, as is the conversion of non-physical [i.e. unperturbed AA^"^) modes 
to real particles when the nuclear density vanishes. A simplified one-dimensional 
scenario is used to investigate the reflection and transmission of pionic modes at the 
nuclear surface. It is found that reflection of pionic modes is rather unlikely, but the 
process can be incorporated into transport descriptions by the use of approximate 
local transmission coefficients. 
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1 Introduction 



In collisions between two heavy nuclei at bombarding energies from a few hundred 
MeV up to several GeV per nucleon, hadronic matter at high density and temperature 
is formed. In such collisions a large number of energetic particles are produced and 
may be used as probes of the hot and dense phase of the reaction |1|, ^. 

Microscopic transport models, such as BUU and QMD |^ ^, ^ |^ have been fairly 
successful in describing particle production in heavy-ion reactions. In these trans- 
port models, the nucleons propagate in an effective one-body field while subject to 
direct two-body collisions. Sufficiently energetic nucleon-nucleon collisions may agi- 
tate one or both of the colliding nucleons to a nucleon resonance, especially A(1232), 
A^*(1440), and A^*(1535). Such resonances propagate in their own mean field and may 
collide with nucleons or other nucleon resonances as well. Furthermore, the nucleon 
resonances may decay by meson emission and these decay processes constitute the 
main mechanisms for the production of energetic mesons 0. 

The transport descriptions normally employ the vacuum properties of the res- 
onances and mesons, i.e. the needed cross sections, decay widths, and dispersion 
relations are taken according to their values in vacuum However, in infinite nu- 
clear matter, a system of interacting vr mesons, nucleons, and A isobars will couple 
to form spin-isospin modes. Some of these modes are non-collective in their charac- 
ter, dominated by a single baryon-hole excitation, while other modes are collective 
and correspond to meson-like states (quasimesons). The non-collective spin-isospin 
modes correspond to those already included in transport simulations by promoting 
a nucleon from below to above the Fermi surface. The collective spin-isospin modes 
can effectively be regarded as particles of mesonic character (quasimesons), which by 
means of a local density and temperature approximation can be incorporated into the 
transport descriptions. 

Some in-medium modifications have already been employed in calculations of 
heavy- ion collisions, both qualitatively and by transport simulations 0, ^, |^. A 
more elaborate vr + NN~^ + AA^~^ model was employed in ref. |10| to derive several 
quantities useful for implementation of in-medium properties in transport descrip- 
tions. While it is straightforward to apply the local density approximation in the 
interior regions of the nuclear system, there are conceptual problems of how to pro- 
ceed at the nuclear surface where the density approaches zero and the quasimesons 
convert to real physical particles. The problem is that while in a real system no hole 
states exist in vacuum, a collective (AA^"^-like) mode can in a stationary infinite 
system exist for an arbitrary small (but finite) density. This is because the particles 
in the stationary system have infinite time to explore the entire system and form 
collective modes also at extremely low densities. 

In this paper we will therefore discuss how a proper conversion from quasimesons 
to real particles at the surface can be performed. Earlier works [|, |^ have treated 
this conversion by various approximations (see further the discussion in section B). 



In this work we will present a somewhat different approach, based on ref. |T^, and 
discuss how the approximations in the earlier works can be improved upon. 

The present paper constitutes a qualitative investigation of treatment of pionic 
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modes at a nuclear surface in transport simulations. To make the presentation simple 
and transparent we will therefore restrict ourselves to some special cases of the more 
complete investigation in ref. |T^. We will only consider the collective modes in the 
spin-longitudinal channel (pion like), since this is the dominant channel at the energies 
we have in mind in this paper. Collective modes in the spin-transverse channel (p 
meson like) can be treated completely analogously. Furthermore we will consider the 
zero-temperature case, and the A width will not be included in the calculation of the 
dispersion relations (denoted as the reference case in ref. ||10|| ). 

The justification for considering only T = in this paper is that there is not 
a strong dependence on the dispersion relations of the collective modes for moder- 
ately large temperatures, especially not at low densities at the surface. Also, the 
temperature is not expected to be very high at the surface. However we want to 
emphasize that there is no principal difficulty associated with incorporating T > in 
the treatment. 

The motivation for omitting the A width, in the calculation of the dispersion 
relations, is somewhat more involved. Including the A width self-consistently in the 
calculation of the spin-isospin modes encompasses decay processes hke 

AiV-i ^ (A^ + ^k)N~^ . 

However, since such processes are already explicitly contained in the transport simu- 
lation by processes like 

+ N ^ N + nk , 

it does not seem to be correct to include the entire self-consistent A width when 
calculating the collective modes to be used in the transport description. Instead it 
seems more correct to use the results obtained with the A width omitted, both for 
the energies of the modes and for the partial A widths to be used in the decays 

A N + TTj. 

However one should note that by omitting the A width in the dispersion relations 
one also fails to take into account the fact that the pionic modes have a Breit-Wigner 
like energy distribution, analogous with the A. The width and center of this distribu- 
tion are determined by the the self-consistent A width and depends on the particular 
pionic mode and its momentum. The center of the distribution approximately corre- 



sponds to the energy found when the A width is omitted |]10 . 

In section ^we will give a brief presentation of the model. The dispersion relations 
and amplitudes of the spin-isospin modes obtained in infinite nuclear matter, are 
presented and discussed in section ^. Section ^ is devoted to a discussion of the 
pionic modes at a nuclear surface and the implications for transport descriptions, 
while our results are summarized in section ^. In addition we present in appendix 
a discussion of reflection and transmission properties at the nuclear surface, and in 
appendix some technical details for the RPA equations. 
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2 The model 



The model presented in this section is treated and motivated in detail in ref. ||T0 |. 
For convenience we here present a brief recapitulation of the essential points. Fur- 
thermore, the presentation in this section only treats the spin-longitudinal channel 
for the special case when T = and the A width is omitted in the calculation of the 
spin-isospin modes. 



2.1 Spin-isospin modes in an infinite system 

We consider a system of interacting nucleons (A^), delta isobars (A) and pi mesons 
(vr). In order to investigate the in-medium properties of the interacting particles, we 
employ a cubic box with side length L; the calculated properties are not sensitive to 
the actual size, so we need not take the limit L — > cxd explicitly. 

The in-medium properties are obtained by using the Green's function technique, 
starting from non-interacting hadrons. The non-interacting Hamiltonian can be writ- 
ten 

Ho = Yl ^kb\K + ^^AQi)^}^i ■ (1) 

k I 

Here the index k = (p^,; Sk,ms^; tk,mtj represents the baryon momentum, spin, and 
isospin. The spin and isospin quantum numbers, Sk and tfc, take the values | and | 
for and A, respectively. The energy of baryon k moving in a (spatially constant) 
potential is denoted e^. The baryon creation and annihilation operators, hi and 6^, 
are normalized such that they satisfy the usual anti-commutation relation, 

{bl k,} = 6,,,, . (2) 

In the pion part of Hq, the index / represents the pion momentum and isospin, / = 
(Pi, A; = 0, ± 1). The meson energy is given by hLU,, = [ml + q'^Y^'^ and the creation 
and annihilation operators of the pion are normalized such that they satisfy the usual 
commutation relation, 

1^1, = ■ (3) 

Note that the A isobar described by Hq has no decay width. Fa = 0. When the 
interactions are turned on, the A width will emerge and it will then automatically 
include also the free width. 



2.1.1 Basic interactions 

At the NnN and A^vrA vertices we will use effective p-wave interactions, Vn-kN and 



V/VttA, which in the momentum representation can be written as [|r^, ^ 



NttN 



IC 



{he 



NttA 



IC 



L3 

ihc) 



2m 



L3 



m^c^ + 
2mA 



In 



NN 



m^c^ 



(4) 



Ina 



F^{q){S+-q,J T • 0^ (g) + h.c. (5) 
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In these expressions, y/s is the center-of-mass energy in the A^vr system and q^^ is 
the pion momentum in the Nn center-of-mass system, which in the non-relativistic 
hmit is given by 

1cm ~ 2 I 9 2 I Pn ' (6) 

mAfC^ + huj rrLNC^ + nuj 
where hu and q is the pion energy and momentum, and pjy is the nucleon momentum 
in an arbitrary frame. The Pauh spin and isospin matrices are denoted cr and r, and 
and T"*" are spin and isospin ^ to | transition operators normahzed such that 
< |, f IS'i |i, i >= 1.0 The momentum representation of the pion field is taken as 



2' 2l +ll2' 2 

<Pl{q) 



+ . (7) 



'2huj.„{q) 

The interactions contain a monopole form factor. 



A2 - {cqf 



and the couphng constants are determined at (cg)^ = {KujY — {cqY — ("^ttC^)^ and 
= rriNC^ or y/s = nij^c^. 

In addition we will also include effective short-range interactions at nucleon-hole 
vertices, again written in momentum space, 

and the corresponding interactions obtained when one (or two) of the nucleons is 
replaced by a A. The strength of the short-range interactions is determined by the 
correlation parameters g'l^j^, o'na^ ^^nd g'^^. 



2.2 RPA approximation 

We want to calculate a spin-isospin mode Green's function within the RPA approxi- 
mation, symbolically 

P; cu) = Go(«, P;cu) + J2 Goia, 7; to) V(7, cu) ^^^^(/t, cu) . (10) 

The spin-isospin modes, here represented by the Green's function G^^^, will in this 
approximation be obtained as an infinite iteration of (non-interacting) pion, nucleon- 
hole, and A- hole states, represented by the diagonal Green's function Gq, coupled with 
the interactions specified in eqs. (§-^ which here are summarized by the symbolic 
interaction V. 

^For clarity, we generally employ bold- face characters to denote quantities with vector and tensor 
properties under ordinary spatial rotations, while arrows are employed to indicate the transformation 
properties under rotations in isospace. 
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In nuclear collisions at beam energies up to about one GeV per nucleon, which 
is the domain of application that we have in mind, only relatively few mesons and 
isobars are produced and so the associated quantum-statistical effects may be ignored. 
Accordingly, we assume tt-a -C 1, and <C 1. 

A set of RPA equations, equivalent to eq. (|10|) were derived in ref. |]T0[ . From these 



equations eigenvectors and eigenenergies are obtained for the different spin-isospin 
modes. The eigenvectors will yield the amplitudes of the different components (vr, 
NN~^, AA^~^) forming the particular spin-isospin mode with the given eigenenergy. 
These RPA amplitudes contain important information about the nature of the differ- 
ent spin-isospin modes. The spin-isospin modes (or excited RPA states), I'if^ >, are 
created by an operator Ql, 

Qliq, A) = E XJkiq, ^)b]k + E Z^Q^ A) vri - 5: W^q, A) rr, . (11) 

jk k k 

The quantity Xj'^(q, A) is here the amplitude of the baryon-hole [NN^^ or AA^~^) 
component of the spin-isospin mode |\E'i, > at momentum q and isospin A, while 
Z^{q, A) and Wj^lq, A), in the same way, are the amplitudes of the pionic component. 
The summation over baryon and meson states in eq. (|lT]) is restricted by taking 
Xjk oc 5p^ p^+q, Zk oc 5p^^q5xk,x, and Wk oc 5p^,^q5x^ -x- 
The RPA equations are obtained from the relation 

< [SQ, [H, gt]] >= hu < [5Q, gt] > , (12) 

with SQ = blbj, ifr, or n}., and where the brackets < ■ > denote the expectation 
value in the interacting ground state. It can be shown that the set of RPA solutions 
constitutes an orthonormal set. For convenience the solutions of the RPA equations 



of ref. |Tn], are recapitulated in appendix 



2.2.1 The total A width 

The A self energy Sa is calculated according to the diagrams in fig. |1], by taking into 
account all the diagrams corresponding to the A decaying into a spin-isospin mode 
and a nucleon, which then again form a A. In the spin-longitudinal channel we obtain 
(ref. m) 

ri(EA,PA) = Im^(^) Y.[9{£)-nip^-q)]M{AN,NA) . (13) 

where the energy available for the spin-isospin mode is given by 

£ = EA-eN{p^-q) , (14) 
and M(34, 12) can be expressed as 



^ | M31;.)M24;.) _ ^ JMk^,;^,,,,, . (15) 
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The factor h{jk, v) is obtained from the interactions at the vertex consisting of 
baryons j and /c, and the spin-isospin mode z/. The interactions to be used depend on 
the non-interacting states that the mode consists of and must therefore be multiphed 
by the amphtude of the corresponding state, 

h{3k,v)'d\3k) = ^^[Z{U) + W{u)]+Y.V,k,mnXmn{^) , (16) 

where d\jk) is a short hand notation for the spin-isospin matrix elements in the spin- 
longitudinal channel, V^-^fc is defined in (^) and (|), Vjk,mn is defined in (|]), and the 
amplitudes X^^^, Z'-, are defined in eq. (PI) . The exphcit expression for h^jk,^) 
is somewhat lengthy and has therefore also been relegated to appendix |B[ 

2.2.2 Specific A channels 

The total A width gives the transition probability per unit time for the A resonance to 
decay to any of its decay channels. In a transport description one explicitly allows the 
A resonance to decay into specific final particles. Consequently, one needs not only 
the total A width (which is the sum of all decay channels) but also the partial widths 
governing the decay into specific RPA channels. These decay channels consist of a 
nucleon and one of the spin-isospin modes. Since we have access to all the amplitudes 
of a given spin-isospin mode on the different unperturbed states, it is possible to 
derive an expression for the partial contribution to from the A decay to a specific 
mode u. The right-hand side of fig. |1| shows a diagrammatic representation of such 
a process. The partial A width for a A decay to a nucleon and a spin-longitudinal 
mode z/ becomes UTOi 



xfiN{p^) (27r)^5(p^ -Pn - q)2'!T6{E^ - Cat - huj) 
= \ j ^ |/i(AiV,u;,)|2n^(p^ - q) 2t^5{E^ - - Uuj,) (17) 

where the factor, = 1 — njv, takes into account the Pauli blocking of the nucleon. 
Note that when this expression is to be used in transport models the factor un should 
be omitted since the Pauli blocking is treated explicitly in the transport description. 
The expression (|T7p is identical to the contribution from one of the u terms in eq. 

m. 



3 Dispersion relations and amplitudes 

From eq. (|59| ) in appendix |B| we calculate the energies of the spin-isospin modes that 
are formed in the interacting system, i.e. their dispersion relations. Fig. |^ displays 
the dispersion relations at normal nuclear density, p = = 0.153 fm~^. In fig. |^ 
a number of different modes in the spin-longitudinal (vr-like) channel are apparent. 
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Some of those are non-collective A^A^ ^ modes (solid curves), which have their energies 
within the regions 

2 

< huj < ^— + — , q<2pF , 

< hu; < ^ + ^, ,>2p^. (18) 

Since we are presenting our results for a box normalization with a finite side length L, 
we obtain a discrete number of non-collective NN~^ modes. The total number of spin- 
isospin modes within the region ([l8| ) depends on L and tends towards a continuum 
in the limit L — > oo. Similarly, a number of non-collective AA^~^ states emerge in fig. 
1^ which, for a fixed q, have their energies constrained to a band, 

m/^— 771^-1 < nuj < — rrij^ -\ 1 . (19) 

2mA T^A 2mA t^a 

The non-collective baryon-hole modes correspond in a transport description to prop- 
agation of uncoupled baryons (A^ or A). This was discussed and studied in detail in 
ref. [|T0| and will therefore not be further discussed in this paper. 

In addition, two collective modes appear in fig. ||, represented by dot-dashed 
curves. The lower one starts at hu = m^c^ at g = and continues into the AA^~^ 
region around q ~ 360 MeV/c. This mode will in the following be referred to as vfi. 
The upper collective mode starts slightly above hu ^ m/\(? — tjinc^ at g = and 
approaches TioJt, = [(m^c^)^ -|- (cg)^]^/^ at large q. This mode is denoted 7f2. 

The incorporation of the two collective spin-isospin modes into transport equa- 
tions is more involved. These modes can be regarded as separate particles of pionic 
character, tti and 7:2, and treated in a manner analogous to the standard treatment of 
the pion. Since the pion is then fully included in the description, it should no longer 
be treated explicitly. The propagation of the two collective pionic modes is governed 
by the effective Hamiltonians 

Hi{r,q) = huji{q; p{r)) = huji , 

H2{r,q) = hu2{q;p{r)) = hOj2 , (20) 

where huji and huj2 are the energy-momentum relations for the lower and upper 
collective modes displayed in fig. ^ for p = Note that the spatial dependence of 
Hi{r, q) is incorporated by representing p{r) as a local quantity. To facilitate center- 
of-mass transformations we parametrize the dispersion relations of the pionic modes 
in the form 

hCo{q; p)^{[cq- cgo(p)]' + mo(p) V}^ + Uo{p) (21) 

For convenience in the transport simulations, we have chosen to use relatively simple 
expressions for the parametrization, rather than to try to optimize the fit. In this 
way a quasipion moves like a relativistic particle with the group velocity determined 
by an effective momentum and energy 

^ = c^^i^ = c^. (22) 
dq huo — Uq huj* 
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The density-dependent parameters qo, mo and Uq are presented in fig. ^. 

Furthermore, in the colhsion term of the standard transport description the pro- 
cess for the production and absorption of pions A ^ + vr, should be replaced by 
the two distinct processes 

A ^ + TTi and A ^ + tt 2 • (23) 

The A decay is governed by the A decay width in the medium to these two specific 
channels, F-^. These partial widths, should be employed in the same manner as the 
free width, i.e. they describe the probability for the A isobar to decay into a nucleon 
and a pion. The only difference is that several collective pionic modes are available 
in the final state. In fig. ^ we present total and partial A widths for a A with the 
momentum 300 MeV/c. The reverse processes in (^) are characterized by cross 
sections that were presented and discussed in ref . . To obtain the partial A width 
from eq. (0) it is necessary to know the amplitudes Z, Z^Xaat-i and J^X^^-i. We 
therefore also present a parametrization of these quantities. On the lower pionic mode 
the pionic component dominates for small momenta, q, and the AA^~^ component 
dominates at larger momenta. Therefore the sum of all individual AA^~^ components 
will for small q increase with q. However, when the lower pionic mode enters the 
AN~^ region the collectivity disappears gradually, and the sum of all individual 
AA^~^ components starts to decrease with q. Thus we employ the form 

$:xlv-i(9;p) ^ /A(g, qxxv-i]) , (24) 

^Xl\.,{q;p) ^ /^(g, C7[X;Vi]) , (25) 

with 

The density-dependent coefficients C are presented in fig. |. On the upper pionic 
mode we instead parametrize the amplitudes as 



5:xxVi(9;p) ^ C'^[^IVi]V/2(?>C''[^IVi]) (27) 



Ex;Vi(g;P) ^ C'^[X^^\.,yi-Uq,C'[Xt,\..]) , (28) 



where 



/2(g,C")= l + exp(C', + C[cq + CU/cq] , (29) 



with the density-dependent coefficients C, displayed in fig. |. The amplitudes Z of 
the pion component are obtained from the parametrization of the squared amplitudes, 
eqs. (gup and (|T]) below, as Z = v^. 

The total A decay width, has apart from the partial contributions F^, also the 
partial contributions F^^ and F^^ . The partial width F^^ gives the probabil- 
ity for the A to decay into a nucleon and a NN~^ state. In a transport description, 
this implies that we initially have a A and after the decay process we have two nu- 
cleons above the Fermi surface and a hole left in the Fermi sea. But this is the same 
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process as if the A would collide with a nucleoli below the Fermi surface to give two 
nucleons above the Fermi surface. This process is normally already included in the 
collision term in a standard transport description, and the probability for such a col- 
lision is given by the cross section for the process A + A^— i>A^ + A^. In a transport 
description it is therefore not correct to both include a A decay according to 
and a collision term with A + — * + A^. Instead, the correct procedure should 
be to exclude F^^ and modify the cross section a{A + N —>■ N + N) to be the 
in-medium cross section. Calculations of such in-medium cross sections was discussed 
in ref. [0. In the same way, F^^ should be excluded in a transport description, 
and ct(A + A^ ^ A + A^) be the in-medium cross section. 

Although the collective pionic modes can thus be effectively treated as ordinary 
particles, the fact that their wave functions contain components from vr, NN~^ and 
AA^~^ states makes it difficult to picture them in a physically simple manner. Fortu- 
nately, their specific structure is not important for the transport process, as as long as 
these quasiparticles remain well inside the nuclear medium. First when such a quasi- 
particle penetrates a nuclear surface and emerges as a free particle is it physically 
meaningful to determine what kind of real particle it is. The gradual transformation 
of the collective quasiparticle occurs automatically within the formalism, because as 
the density is lowered, p — 0, a pionic mode will acquire 100% of either the pion com- 
ponent or the AA^~^ component, depending on uj and q. That is to say, it will turn 
into either a free pion or an unperturbed AA^^^ state. The squared amplitudes have 
their values between zero and unity and we therefore employ the parametrizations 

^*H9;p)' ^ /2(g, (30) 

Z^^{q-pf ^ l-h{q,C"[Z^^]) (31) 
^ l-/2(g,C7"[XlVi]) (32) 
Y^Xl^^.,{q;pf ^ /2(g, (7"[XXVi]) , (33) 

with f2iq,C'') from eq. (^), while the sum of the squared NN~^ amplitudes are 
obtained from the normalization, i.e. all squared amplitudes sum up to unity. The 
density-dependent coefficients C" , are presented in fig. 0. There remains the practical 
problem of how to represent an unperturbed AA^~^ state when p — >■ 0. However, as 
will be discussed in section ^, only a very small fraction of the pionic modes (vanishing 
for a stationary density profile) will emerge as unperturbed AA^~^ states. 

Note that this approach is different from earlier works 1^, ^ where the nature 
of the pionic mode was determined already in the creation process, i.e. when the mode 
was created it was determined whether it represented a free pion or an unperturbed 
AA^^^ state. This difference in approach will have crucial effects for the collective 
modes that escape the system, as will be seen in the next sections. 



4 Quasipions at the nuclear surface 

In previous works the quasipions at a surface have been treated in various approx- 
imations. In ref. effective dispersion relations were introduced, corresponding to 
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modes with either 100% pionic or AA^~^ component. The pionic mode was then prop- 
agated as a quasipion, emerging as a free pion at the surface, while the AA^~^ mode 
was used to derive a A potential for the uncoupled As. In this way only pions and 
As escape the system, but the drawback is that the effective dispersion relations are 
quite distorted compared to the original ones, and that the in-medium effects seem 
to be over-estimated by allowing all uncoupled As propagate in a collective potential. 

In both collective modes were propagated, and hence some modes escape the 
system as unperturbed AA^~^ states. This was effectively taken care of by converting 
these modes to free As (neglecting baryon number conservation), with the justification 
that the number of modes escaping the system as unperturbed AA^~^ states were 
found to be small. We will in this paper report on alternative ways to treat the 



modes at the surface based on ref. |]T0[, and how this treatment can improve the 
descriptions in refs. [|], 

We first consider the simplified case when the nuclear surface is stationary i.e. 
p{r, t) = p{r, 0) for all times t. This means that there is no explicit time dependence 
in the effective Hamiltonians in equation (ppj), and thus the energies of the collective 
pionic modes are conserved. For a pionic mode with energy hujy{q\ p) propagating in 
a varying density this means that the momentum q will change as p changes. That 
is to say, the pionic mode effectively feels a potential. In addition to momentum 
changes, also the amphtudes Xji^{q; p), Z'^{q; p) and W^{q; p) of the baryon- hole and 
pion components will change as the density changes. In the limit of vanishing density 
either X'^^-i or Z'^ will turn to unity, depending on the mode v and its energy and 
momentum. 

In the latter case no problems arise, the collective mode has simply been con- 
verted to a free pion escaping the system. However, in the former case there is some 
inconsistency in the formalism since there are no hole states in vacuum (p = 0). In a 
quantal description of spin-isospin modes propagating at a surface different scenarios 
could emerge. There is some small probability that the mode could be reflected at 
the surface. Alternatively the mode could break up in an uncoupled A escaping the 
system, with the hole is trapped inside the nucleus. In a transport description this 
could be handled by allowing the test particle representing the mode v to absorb a 
nearby nucleon, converting it into a A isobar. 

We consider in this section the quasipions at a surface propagating from normal 
nuclear density to vacuum. In this case the amplitude will be very small and we 
will therefore omit it in the qualitative discussion of this section. 

In fig. 1^ we present the energy Tiuji {q) of the lower collective mode tti for different 
densities in the range 0.1p° < p < p°. The line closest to the free pion relation 
(dotted line) represents the dispersion relation at the lowest density. As the density 
is increased the energy relation is (for each fixed q) lowered. For small q values the 
mode is completely dominated by the pion component. As q is increased also the 
AA^^^ components starts to contribute substantially (although to different extent 
depending on the density). In the approximate range 300 MeV/c < q < 400 MeV/c 
(depending on density) the mode enters the AA^~^ continuum and changes character 
from collective mode to non-collective. 

We will in this section discuss four different examples, suitably chosen to illustrate 
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the main features. 



4.1 Lower pionic mode 

In our first example we consider the mode tti created at normal density with energy 
huji = 200 MeV propagating towards vacuum without any interactions. Initially this 
mode has the following characteristics: 

g^ 220 MeV/c \Z\^ ^ 0.72 Eatv-i I^Aiv-i P ~ 0.25 Etvjv-i l^iv^-^ P ~ 0.03 . 

As the density decreases it will, due to the energy conservation, follow the path 
indicated by the dashed line in fig. ^a, from q ^ 220 MeV/c to g ~ 150 MeV/c. In 
fig. we see how the squared amplitudes vary with the density for this particular 
energy of the mode tti. As seen in fig. only the pion component remains, as the 
zero density limit is reached. Thus in this particular case, the mode will escape the 
system as a free pion. In fig. p!0| a we show the squared amplitudes of the individual 
AA^~^ amplitudes for the energy hui = 200 MeV of the mode tti. We see that the 
mode is collective at all densities. 

In the discussion so far we have not addressed the creation process. How probable 
is it that we produce the mode tti at the particular energy hui = 200 MeV? This 



information is given by the partial A decay width to the mode tti. In fig. 11a we 



display the partial A decay width, F^, for different densities, as a function of the 



energy of the emitted quasipion vfi. Note that the width displayed in fig. |Tl|a is for a 
A at rest, and that the Pauli blocking of the emitted nucleon has not been taken into 
account (the Pauli blocking is treated explicitly in a transport description). As seen 
in fig. pA| a the width is quite substantial at the quasipion energy hoji = 200 MeV, 
about 80 MeV compared to the free width, which is about 30 MeV at this energy. 
Thus it is quite probable for a A to decay to the mode tti around this quasipion 
energy. 

Apart from penetrating the surface there is also some probability for the mode tti 
to be reflected at the surface. This is difficult to exactly predict since the reflection 
coefficient, TZ, will depend on the actual density profile. However a first estimate 
can be obtained by considering a one-dimensional scenario with a density profile 
that corresponds to a Wood-Saxon potential with a surface thickness a = 0.65 fm. 
The height of the potential is then given by the change in the momentum of the 
pionic mode, and for the case of hui = 200 MeV we obtain TZ ~ 0.017, see further 
the discussion in appendix ^ This number is very small, which implies that the 
reflection can practically be neglected for this particular case. However, in other 
situations where the momentum change is different or the density profile is sharper, 
the reflection coefficient may be larger. We have therefore devoted appendix ^ to 
discuss how effective local reflection coefficients can be obtained and implemented in 
transport descriptions. 

In our second example we also consider the mode vfi, but now with energy huji = 
295 MeV. Initially this mode will now have the characteristics: 

g^ 450 MeV/c \Z\^^0.0 Eatv-i I^atv-i P ~ 0.94 Ejvtv-i l^^iv-i P ~ 0.06 . 
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As the density decreases it follows the path indicated by the dot-dashed line in fig. ^a, 
from q k, 450 MeV/c to g ~ 250 MeV/c. In fig. we see how the squared amplitudes 
vary with the density for this particular energy of the mode tti. As seen in fig. ^ also 
in this case only the pion component remains, as the zero density limit is reached, 
although TTi initially was completely dominated by the AA^~^ component. Thus also 
in this case the mode will escape the system as a free pion. 

In fig. pi]| & we show the squared amplitudes of the individual AN~'^ amplitudes 
for the energy huji = 295 MeV of the mode tti. Here the situation is different from 
the case at huji = 200 MeV, since initially the mode is dominated by a single AN~^ 
component, i.e. the mode is non-collective. However as the density is lowered the 
strength is spread over more AA^^^ components and at low densities the mode is 
completely collective. 

From fig. |ll]a we see that the partial width is very close to zero at the quasipion 
energy hoji = 295 MeV (dot-dashed curve) at normal nuclear density. Thus a A at 
normal nuclear density will not decay to a quasipion with this energy. The partial 
A width becomes very small because the collective strength, as well as the pion 
component, is negligible in this case. Although the mode at this energy cannot be 
created at normal nuclear density, it may still be created at lower densities, as can 
be seen in fig. pl]a. 

Making the same assumptions as in the first example, we find that the refiection 
coefficient, becomes smaller than 10~^ for this case. 

We thus conclude that the lower pionic mode penetrating the surface will always 
emerge as a free pion. For low energies this is natural because as the density decreases 
the dispersion relation of the pionic mode approaches the free pion relation. The cases 
of sufficiently large energy to approach the unperturbed AA^~^ branch as the density 
is lowered will never occur since no such modes will be created from a decaying A, 
because the partial A width will be zero. 

4.2 Upper pionic mode 

Also on the upper collective mode, 7^2 a similar effect will occur. However some 
properties are somewhat different so we will therefore illustrate also this case with 
two typical examples. In our third example we thus consider the mode 712, with energy 
hLJ2 = 320 MeV. Note that at this energy the mode can only exist at densities up to 
about 0.5p°. Initially, at p = 0.5p'', this mode will have the characteristics: 

q^OMeV/c \Z\^^0.0 Ean-^ I^an-A^ ^ l-O Enn-^ I^nn-A^ ^ 0.0 . 

As the density decreases we follow the path indicated by the dashed line in fig. 
from g ^ MeV/c to g ~ 270 MeV/c. Note that on the upper collective mode the 
momentum increases as the density decreases, corresponding to a negative potential 
step. 

In fig. ^ we see how the squared amplitudes vary with the density for this par- 
ticular energy of the mode 7i2- As seen in fig. ^c, contrary to previous examples, 
only the AA^~^ component remains, as the zero density limit is reached. In fig. pUj c 
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we show the squared amphtudes of the individual AA^ ^ amphtudes for the energy 
tiuj2 = 320 MeV of the mode n2- Initially the mode is collective, but as the density is 



lowered the mode becomes more and more non-collective. From fig. |TT^ we however 
see that the partial width T\ actually is zero at the quasipion energy hu2 = 320 MeV 
(dot-dashed curve) at all densities. Thus a A will not decay to a quasipion with this 
energy. 

A mode at this energy with a very low momentum at half nuclear density, which 
could occur in a time- dependent density, could have a very large reflection coefficient, 
approaching unity as the initial quasipion momentum approaches zero. 

In our fourth and last example in this section we consider the mode 7f2 with energy 
huj2 = 380 MeV. Initially this mode will now have the properties: 

g^l70MeV/c \Z\^ ^ 0.12 Eatv-i I^Aiv-^ I' ~ 0.88 Etvjv-i |^7Viv-i P ~ 0.0 . 

As the density decreases we follow the path indicated by the dot-dashed line in fig. 
from q ^ 170 MeV/c to g ~ 350 MeV/c. In fig. ^ we see how the squared amplitudes 
vary with the density for this particular energy of the mode 7f2. As seen in fig. ^ 
in this case only the pion component remains, as the zero density limit is reached, 
although #2 initially was dominated by the AA^~^ component. Thus also in this case 
the mode will escape the system as a free pion. 

In fig. [To| (j we show the squared amplitudes of the individual AA^^^ amplitudes 
for the energy huj2 = 380 MeV of the mode 7^2. We see that the mode at all densities 
is completely collective. From fig. we see that the partial width F^ is quite 
substantial at the quasipion energy hLJ2 = 380 MeV (dot-dashed curve) also at normal 
nuclear density. The reflection coefficient, making the same assumptions as in the first 
example, becomes for this case smaller than 10~^. 

We thus conclude that also the upper pionic mode penetrating the surface will 
always emerge as a free pion. For high energies this is natural because as the density 
decreases the dispersion relation of the pionic mode approaches the free pion relation. 
The cases of low energy when the unperturbed AA^~^ branch is approached as the 
density is lowered, will never occur since no such modes will be created from a decaying 
A, because the partial A width will be zero. 



4.3 Refined scenarios 

If the surface changes with time the energy of the pionic mode need not to be con- 
served, and there is some small possibility for the pionic mode to end up as an unper- 
turbed A-hole state. The actual fraction of such modes is hard to estimate without 
an explicit transport simulation, but based on the scenario for the time-independent 
density profile, it is reasonable to expect that only a very small fraction of the pionic 
modes will end up as unperturbed A-hole states in vacuum. 

In a quantum description such a mode could be either reflected at the surface, or 
the mode could break up into an uncoupled A and hole, where the hole is trapped 
inside the nucleus, and the A escape the system as a free A. Based on the results 
discussed in the explicit examples of this section, and the presentation in appendix 
we expect that the reflection at the surface will be very small. 
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In a transport description the reflection can be incorporated by a local transmis- 
sion coeflicient as discussed in appendix 0. If the pionic mode is not reflected at the 
surface, and its amplitude approaches 100% of the AN~^ component as the density 
approaches zero, the mode should thus break up into an uncoupled A and hole. In 
a transport simulation this could be practically handled by allowing the pionic mode 
to absorb a nearby nucleon, forming an uncoupled A, when the density falls below a 
specifled value. This prescription has some quantum mechanical justiflcation by the 
fact that at very low density a wave packet representing the pionic mode, will not be 
very well localized, but instead have a large spatial spread. 

5 Summary 

In-medium properties obtained in an inflnite stationary system consisting of inter- 
acting nucleons, nucleon resonances and mesons, can be incorporated into transport 
descriptions by a local density approximation. While such a prescription is rather 
straightforward to implement in the interior regions of the nuclear system, concep- 
tual problems exist at the nuclear surface. When the nuclear density approaches zero, 
collective mesonic modes formed in the medium have to be converted to real particles 
in vacuum. The problems arise since some collective modes (e.g. AA^^^-like) may 
exists in the inflnite stationary system at arbitrary low (but non-vanishing) density, 
but no corresponding real particle exists in vacuum. This problem has been apparent 
in previous works ||^, Q where collective modes have been incorporated into transport 
descriptions as quasimesons. The character of the quasimesons {i.e. realization in 
vacuum) were in those works determined already at the time of creation. 

Based on the formalism of ref. |]10|, we have in this paper employed a more elabo- 
rate 7r + NN~^ -|- AA^~^ model (relative to the works H, §) to investigate a somewhat 
different treatment of the collective pionic modes at a nuclear surface. In this formal- 
ism we have obtained not only density dependent dispersion relations of the pionic 
modes, but also density dependent amplitudes of the components constituting the 
pionic mode. These quantities are conveniently parametrized with density dependent 
parameters, in section ^ 

For the transport process it is not needed to determine the character of the pionic 
modes until they penetrate the surface and emerge as free particles. This is auto- 
matically determined within our formalism from the amplitudes at zero density. We 
have further showed in section |^ that for a stationary density proflle, the conserva- 
tion of the energy of the pionic mode and the partial A decay width, together leads 
to the fact that only real pions are realized as free particles when the pionic mode 
penetrates the surface. Note that this flnding is different from earlier works, and it 
demonstrates the importance of deriving dispersion relations and partial A widths 
consistently within a realistici AA^~^ model. 

In a more reflned scenario, where the density changes with time, deviations from 
this picture can be expected and also the unperturbed AA^~^ component may be 

^By "realistic" we here mean that the AA^^^ model contains a continuum of AA^^^ and NN~^ 
states, as compared to the more simple two- level AA^^^ model used for example in refs. [§[ p). 
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realized in the limit of vanishing density. The actual fraction of such modes is hard 
to estimate without an explicit transport simulation, but based on the arguments in 
section ^ for the stationary surface, we expect this fraction to be very small. 

For the rare cases when the unperturbed AA^~^ component is realized in the 
limit of vanishing density, the pionic mode must be converted to real particles. In 
an extended description this could be made by allowing the mode to break up into 
an uncoupled A and a A^~^, where the hole remains trapped in the nuclear system 
and the A escapes the system. Based on our formalism, this seems to be the most 
probable scenario, since the collective strength disappears on these modes as the 
density approaches zero. This could be implemented in a transport simulation by 
letting the pionic mode absorb a nearby nucleon to form an uncoupled A, as may be 
justified by the quantum-mechanical feature that the wave packet representing the 
pionic mode is not well localized at very low density. 

Alternatively the pionic modes could in a quantal description be reflected at the 
surface. We have investigated the reflection and transmission probabilities for the 
collective modes in a simplified one-dimensional scenario, where the modes propagate 
perpendicular to the surface. Exploring typical scenarios, we have found that (with 
only very few exceptions) the reflection of the pionic modes will be smaller than a 
few percent, however, in appendix ^ we have suggested how the reflection and trans- 
mission probabilities could be incorporated into the transport descriptions by using 
approximative local transmission coefficients. 

Stimulating discussions with Volker Koch are acknowledged. This work was supported 
by the Swedish Natural Science Research Council, and by the Director, Office of En- 
ergy Research, Office of High Energy and Nuclear Physics, Nuclear Physics Division 
of the U.S. Department of Energy under Contract No. DE-AC03-76SF00098. 

A Reflection and transmission at a surface 

In the semi-classical transport descriptions all (test)particles are treated as classical 
particles. When such a classical particle propagates in a spatially varying potential 
its velocity will change. But assuming its energy exceeds the maximum value of the 
potential, the particle will continue to propagate through the varying potential. This 
is in contrast to a quantal description where the particle is represented by a wave 
packet. This wave packet has some probability to be reflected in a spatially varying 
potential. The reflection probability depends on several factors, such as the energy 
of the wave packet, and the height and shape of the potential. 

In this section we investigate reflection and transmission probabilities in differ- 
ent idealized situations for a one-dimensional scenario, corresponding to the direction 
normal to the nuclear surface. We will also discuss how these effects could be ap- 
proximately incorporated into transport descriptions. The treatment is intended to 
be used for the collective spin-isospin modes (quasi-pions) approaching a nuclear sur- 
face. The formalism, though, is quite general and could be applied also for other 
particles. 
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We will start to investigate some special cases when the potential is stationary. 
Subsequently we will discuss how the non- stationary case could be treated. In all 
cases we will assume that the potential is constant outside a finite interval, 



X < xl 
V{x) = { V{x) Xl < X < xr 

Vo Xr< X 



(34) 



For X < a;^ we have an incoming wave and we ask for the probability that we have 
an outgoing wave at a; > xr, i.e. we seek the transmission and refiection coefficients, 
T and TZ respectively, with T + TZ = 1. Only for some special simple forms of the 
potential V{x) can the coefficients TZ and T be obtained analytically. 
Note that both the stationary Schrodinger equation 



2m dx'^ 



+ V{x) 



ip{x) = Etp{x) 



(35) 



and the stationary Klein-Gordon equation 



i-2 2 " I i-2 2 2 4 

h c -— + Ti uj — m c 
dx^ 



Tpi^x) = 2hujV{x)ip{x) 



(36) 



for a particle with energy E = hu, can be rewritten in the form 

V^(x) = . 
by introducing a local wave number 



or 



k{x) 



k{x) = ^2m[E -V{x)]/h , 



rcj2 - m'^c^ - 2hu;V{x)]/hc 



(37) 



(3^ 



(39) 



respectively. The results obtained in the sequel in this section are derived using the 
Schrodinger equation, but from eqs. (^) to (|37| ) it follows that the results are valid 
also for a particle described by the Klein-Gordon equation, if the local wave number 
is taken according to equation ( |5D| ) instead of equation (P8|). 
The simplest case is a potential step at x = 0, 



V{x) 



X < 
Vo < X 



(40) 



This case can be found in elementary textbooks on quantum mechanics. Here we 
briefiy recapitulate the main steps as a preparation for more complicated cases. The 
Schrodinger equation has the solution 



ip{x) 



_A_^ikLX _j_ 
C ^iknx _|_ 



B -ikr.x 
D_p-ikRX 



X < 
< X 



(41) 
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where ki = V2mE/h and kn = y2m(£' — Vo)/h, for a particle with energy E > Vq. 
Taking only an outgoing solution at x > xr {i.e. D = 0), and relating the coefficients 
A, B and C by a smooth joining at a; = 0, 

r ^(0 -e)= ^(0 + e) 
the reflection and transmission coefficients are obtained as 



n 



W {kL + kuf 



T = ^ = (43) 

Also for other special forms of the potential can the reflection and transmission 
coefficients be derived analytically, such as for the Woods-Saxon type of potential, 

r(x) = • (44) 

Here the reflection coefficient TZ is given by |T3[ 



and the transmission coefficient is obtained from T = 1 — TZ. Note that the case of 
the potential step, eq. (plD, emerges in the limit a — 0. 

For an arbitrary potential the coefficients TZ and T do not have an analytical form. 
An approximate solution can be obtained by approximating the potential V{x) by a 
piecewise constant potential 

{0 X < Xq = Xl 

V{x)^Trj=lVj{x) Xo = XL<X<Xn = XR , (46) 
Vo Xn = XL < X 

with 

V,{x) = V{x,+Ax/2Mx,^,-x)-e{x,-x)l I ^^l^i^y^ , (47) 

and by making the ansatz, 

n-l 

i^(yX) ~ ^ 4'j{x)[9{Xj^i — x) — 9{Xj — x)] Xq = Xl < X < Xn = Xr , (48) 

i=o 

with 

^.(x) = -^e'^^'^ + ^e"*^^" (49) 
Jkj Jkj 
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and 

kj = ^2m[E - V{xj + . (50) 

The coefficients Aj and Bj are determined from the condition of smooth joining, 

/ V'j(a;j+i - e) = V^j+i(xj+i + e) ^ _^ n+ 

and the boundary condition of an outgoing solution at z > zr, i.e. An-i = 1 and 
Bn-i = 0. Not that in the especially simple case of no reflection, the ansatz above 
becomes identical to the WKB approximation, in which the the coefficients Aj and 
Bj are constants, i.e. Aj = A and Bj = B for all j. 
The transmission coefficient, T, is obtained from 

^^w- '''' 

and TZ = 1 — T. We have numerically checked that for the Woods-Saxon type of 
potential above, this approximation procedure, with a very high accuracy, yields 
the same reflection and transmission coefficients as the analytical result. We have 
compared the analytical and numerical results for many different energies (E) and 
parameters of the potential (Vq, a), and only for E very close to Vq are there deviations 
larger than a few percent. 

In a heavy-ion collision the potential is not stationary, but changes with the time. 
The formalism above can be generalized to a time-dependent potential. The particles 
to be propagated are in a quantal description represented by wave packets. We 
therefore generalize the definition of the transmission coefficient to be valid also for a 
propagating wave packet in a non- stationary potential. As before we assume that we 
have an incident wave or wave packet, ipmdx, t), ed, x = x^. We associate a probability 
current with this wave packet by 

Jinc(x, t) = - Im ^^^^^^y] (53) 

m \ ax J 

Similarly we assume that we have an outgoing wave or wave packet, ipo^t{x,t), at 
X = xr, and we associate a probability current, jont{x,t), with this wave packet. 
A transmission coefficient can then be defined as the ratio between the outgoing 
probability flowing through the point xr and the incident probability flowing through 
the point xl, 

II^ jontixR,t)dt 
Jt\' jmcixL,t)dt 

where the times ti and t2 are chosen such that a wave packet will propagate through 
both the points xl and xr within the time interval. Note that this definition of T 
agrees with the results stated for the stationary case. By solving the time-dependent 
Schrodinger equation numerically, and expanding a wave packet in this basis set, a 
transmission coefficient could in principle be obtained from eq. (E^ for an arbitrary 
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(known) potential V{x,t). In a transport description, however, the potential V{x,t) 
is not known in advance, and the described method is not so useful. Instead it 
would be more useful with a local transmission coefficient t{xj), which would give 
the probability for the wave packet to be transmitted from the point xj to the point 
Xj+i = Xj + Ax, where the total transmission coefficient is given by T = 11^=0 ^(^^j)- 
Such a local coefficient can formally be obtained by writing 

r = ^-LJ^(^)'^ ^ (55) 
and by approximating the integral in the exponent by a discrete sum, 

_j n—1 n—1 

T ^ e-S^-o = n e-'-^"^)^" = n tixj) . (56) 

j=o j=0 

However, also in this expression the local reflection coefficient r{x) will depend on 
the particle energy, and the shape and magnitude of the potential, and it is therefore 
not straightforward to implement a local transmission coefficient into a transport 
description. 

It is therefore our strategy to approximate the function r{x) with a relative simple 
expression depending on a few parameters, that are adjusted such that eq. (|^) will 
yield a good approximation to the total transmission coefficient for a large class of 
realistic particle energies and potential parameters. We have empirically found that 
the ansatz 



Ci / din k{x) 



2 



with Ci ~ 0.116, yields a good approximation for a range of different particle energies 
and a range of different values of a and Vq in the Woods-Saxon type of stationary 
potential. In fig. ^ we present some examples of this approximation, compared to 
the analytical solutions. 

As was discussed in section ^ it is only in rather few and special cases that the 
transmission coefficient for the collective pionic modes will deviate substantially from 
unity. We therefore conclude this section by summarizing that the reflection of pionic 
modes can be approximately incorporated into transport descriptions by an approxi- 
mate local transmission coefficient 

The error in this approximation is in most realistic cases smaller than 10%. Those 
events where a larger error may occur for pionic modes at a nuclear surface are ex- 
pected to be very rare in transport simulations. Thus the total error in a transport 
treatment by using a local transmission coefficient according to eqs. (|56|) and (^81) 
should be very small. However, the importance of incorporating the effects of reflec- 
tion of pionic modes at the surface seems to be quite small, since in most events the 
total transmission coefficient will be close to unity. 
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B RPA equations for spin-isospin interaction 



In this section we recapitulate the RPA equations, derived in ref. |]T0[, in the spin 
longitudinal channel for the case when T = and the A width is omitted in the RPA 
equations. 

The spin-isospin excitations are characterized by the momentum q and the isospin 
A. The energy of the spin-isospin modes are obtained from the determinant of the 
equation 

1 ^^^N^N^N 1 _^ -j^AA^y^A^A 



X 



N 







(59) 



while the amplitudes are obtained from the solution of eq. 

1 



Z 

w 



1 



and 



3/2 



x{tj,tk; q, A, u) 



^jkiq, -A) 



(60) 
(61) 

(62) 



with 



x(l/2,l/2; q; A, cu) 
x(l/2,3/2; q, A, u) 



X 



N 



x(3/2,l/2; q; A, cu) 



X 



x(3/2,3/2; q, A, w) = 



and analogously for other quantities. In the above expressions a,l3 = N, A, and we 
have used the notations 



J NaJ Nf3 



m\'9L+R?R^\Fj\cq,rD:], 



tRtF^ 



yl3 



2mQ,c^ 



(63) 
(64) 

(65) 

(66) 
(67) 

The numerical factors M.^ = A and M.^ = 16/9 originates from the spin-isospin 
summation, and we have defined the Lindhard functions 



{En{Pn) + T^^f - {cq + cp^f 
{nij^c? + Re huY — c{qY = ^/sl 



UInC 



<l>^(u;,q) 



E 



n{p) - n[p + q) 



L) {p + qy/2m*j^ - p'^/2m%i -huj 



(69) 
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n[P) ^ Hp) 



(70) 



with 



6e 



{piy p' 



AN 



Se 



2mA 2mJ^ 
iPl? P' 



AN 



2mA 2m 



+ Am + V/(£+,p+)-^c^ 
+ Am + Vf{e^,pl) + hto 



(71) 
(72) 



N 



P^ 



Am = mA - mAr; e± = ttiat + ± /io;; p J = p ± g . (73) 

The quantities and are determined from the normahzation of the RPA 
states 



±1 = M^vnx*nXn + M'^VAxlxA + Z*Z -W*W 



1 + }^NN^N(i,N 
y^NAj^A^A 

1 



(74) 



N„.Nif,N -L -h KV « 



y\>NA 



with 



VN{^^,q) 

VA{i^,q) 



he 



E 



n[p) - + q) 



d 



$iv(c^, q) 



'^1 E 



n{p) 



d 



$A(^,g) 



(75) 
(76) 



and where we have assumed dV^ /dfiu ^ 0. 

The factor h{jk\ u), motivated in eq. ([iBD , can using the RPA solution be exphcitly 
expressed as 



h{jk,u;uj,q) 



u)Xa\y^u ) 

a=N,A 



where 



Xtt fir 

vb i^) q) - 9a,jk q) 



(77) 
(78) 
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niN = 940 MeV/c2 
mA = 1230 MeV/c2 

= 140 MeV/c^ 
Ag = 1.5 GeV 


9'nn = 0.9 
g'NA = 0.38 
9'aa = 0.35 
A'^ = 1.0 GeV 


Inn = 1.0 
fivA — 2.2 
/aa = 
Po = 0.153 fm"^ 


V^~Vn = 25.0p/po MeV 


m^ = m^/[l + 0.4049(p/po)] 



Table 1: Parameter values used in the numerical calculations. 

Figure 1: 

Diagrammatic representations of the A self energy (left-hand side) and the partial 
A decay width to the particular spin-isospin mode u, (right-hand side). 

Figure 2: 

The dispersion relations for the spin-isospin modes in the spin-longitudinal channel, 
in infinite nuclear matter at normal nuclear density and zero temperature. The non- 
collective modes are shown by solid curves, while collective modes vfi and fi2 are 
represented by dot-dashed curves. As a reference, the free pion dispersion relation 
hajTril) = [(^TiirC^)^ + {c^yY^'^, and the unperturbed AA^~^ relation hu^^-i^q) = 
m^c^ — uinc^ + g^/2mA, are included as a dotted curves. 

Figure 3: 

The density- dependent parameters Uq (solid curve), go (dashed curve) and mo (dot- 
dashed) used in eq. ( ^1]) for fitting the dispersion relations of the lower collective 
mode TTi (part a) and the upper collective mode (part h). 

Figure 4: 

The total A width F^* and its partial contributions from different spin-isospin modes 
for a A with momentum 300 MeV/c. The solid curve represents the total width, 
the long-dashed line is the contribution from the non-collective NN~'^ modes, the 
short-dashed line is the contribution from the non-collective AA^~^ modes, the dot- 
dashed line is the contribution from the lower pionic mode tti, and the dot-dot-dashed 
line is the contribution from the upper pionic mode 7f2. The intermediate nucleon in 
the A decay is not Pauli blocked. The error bars indicate the estimated uncertainty 
associated with the classification procedure (see ref. [p!0[] ). 

Figure 5: 

The density-dependent parameters, Ci (dot-dashed curve), C2 (short-dashed), C3 
(dot-dot-dashed) and C4 (long-dashed), for the amplitudes of the AA^~^ (part a) and 
the NN~^ (6) components on the lower pionic branch, eq. (0) and (125|). 

Figure 6: 

The density-dependent parameters, Cq (part a), C[ (6), C'_i (c) and C'^ (d), for the 
amplitudes of the AA^~^ (long-dashed curve) and the NN^^ (dot-dashed) components 
on the upper pionic branch, eqs. ( P7| ) and (P^). 
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Figure 7: 

The density-dependent parameters, Cq (part a), C" (6), and C'!_i (c), for the squared 
amphtudes of the vr (dot-dashed curve) and the AA^~^ (short-dashed) components 
on the lower pionic branch, eqs. (pOD and (|3^ ) and of the vr (dot-dot-dashed) and the 
AA^~^ (long-dashed) components on the upper pionic branch, eqs. ( ^1]) and (|33|) . 



Figure 8: 

The energy huj{q) for different densities in the range O.lp" p < ■ In part a is 
presented the energies of the lower collective mode tti while the energies of the upper 
collective mode 1x2 is shown in part h. As references the free pion dispersion relation 
and the unperturbed AA^~^ relation have been plotted as dotted lines. The dashed 
and dot-dashed horizontal lines indicate the energies considered in the four examples 
of sections 0| and O. 



Figure 9: 

Sum of all squared amplitudes of the AA^~^ and NN~^ components, as well as squared 
amplitude of the pion component. In part a and h is shown the squared amplitudes 
for the lower pionic mode vfi at the energies hCoi = 200 MeV and hui = 295 MeV, 
respectively. In c and d is displayed the squared amplitudes for the upper pionic 
mode TTi at the energies hLJ2 = 320 MeV and hoj2 = 380 MeV, respectively. The dot- 
dashed curve represents the pion component, the AA^^^ component is represented by 
a short-dashed curve, and the long-dashed curve represents the NN~'^ component. 



Figure 10: 

Squared amplitudes of the individual AA^^^, components. In part a and b is shown 
the squared amplitudes for the lower pionic mode tti at the energies huji = 200 MeV 
and hoji = 295 MeV, respectively. In c and d is displayed the squared amplitudes 
for the upper pionic mode 712 at the energies hoj2 = 320 MeV and hLJ2 = 380 MeV, 
respectively. 



Figure 11: 

Partial A decay width for a A at rest, for different densities in the range 0.1p° < 
p < p^. In part a is presented the partial A width for the lower collective mode tti 
while the partial A width for the upper collective mode 7r2 is shown in part b. As a 
reference the free A width has been plotted as a dotted line. The dashed and dot- 
dashed vertical lines indicate the energies considered in the four examples of sections 



4.1 and 4.2 



Figure 12: 

Total transmission coefficient, T, obtained for a Woods-Saxon type of potential, eq. 
(PI), for a range of potential parameters Vq and a. The solid curve represents the 
exact expression, eq. (^5]) , while the dashed curve represents the approximative result 
obtained from the local transmission coefficient, eqs. (EBI) and (B3). 
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